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We study the mechanisms of frequency synchronized cluster formation in coupled non-identical 
oscillators and investigate the impact of presence of a leader on the cluster synchronization. We find 
that the introduction of a leader, node having large parameter mismatch, induces a profound change 
in the cluster pattern as well as in the mechanism of the cluster formation. The emergence of a leader 
generates a transition from the driven to the mixed cluster state. The frequency mismatch turns 
out to be responsible for this transition. Additionally, for a chaotic evolution, the driven mechanism 
stands as a primary mechanism for the cluster formation, whereas for a periodic evolution the 
self-organization mechanism becomes equally responsible. 
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The interaction between individual units of a system 
leads to many emerging behaviors, among which synchro¬ 
nization is one of the most fascinating phenomenon which 
has been gaining tremendous attention since the first ex¬ 
perimental demonstration of the phenomena by Chris¬ 
tian Huygens [l|. The unexpected sway and twist of the 
Millennium bridge, synchronization of the neurons in the 
brain and synchronous fire flies are few examples 

of the synchronization in the real world systems. Syn¬ 
chronization is defined as an emergence of some relation 
between the functional of two processes due to interac¬ 
tion Q. Earlier studies on coupled non-identical oscilla¬ 
tors have shown that the exact synchronization is hard 
to achieve when there is a parameter mismatch in the 
local dynamics rather they exhibit phase or general¬ 
ized synchronization Additionally, there exists a 

nontrivial transition to the global phase synchronization 
in a population of globally coupled chaotic non-identical 
oscillators Q- Moreover, cluster pattern and frequency 
of the nodes in a cluster have been shown to be controlled 
by local external forcing as well as by changing the net¬ 
work architecture Q. Furthermore, the coupled oscilla¬ 
tors with heterogeneous coupling have been reported to 
exhibit synchronization triggered by the oscillators hav¬ 
ing strong couplings, further facilitating the synchroniza¬ 
tion among the nodes having weak coupling [9| . 

We present results of cluster synchronization due to its 
importance and occurrence in various real world systems 
represented in terms of interacting units M- We study 
mechanisms of formation of frequency synchronized clus¬ 
ters in coupled non-identical oscillators and investigate 
the influence of a leader on the dynamical evolution of 
other nodes. One possible way of defining a leader is to 
make the natural frequency of a node much higher than 


*sarika@iiti.ac.in 


that of other nodes in the network [TJ, [T^. This is one 
of the traditional way to define a leader in the coupled 
dynamics on network models. Some other ways of defin¬ 
ing leaders are those which depend upon the application 
and motivation of the problem, e. g., Ref.[l^ considers a 
leader which exchanges information with its neighbors as 
well as has access to its own state. In Ref. the neuron 
which fires first is considered as leader. Further, a leader 
can be assigned based on its degree in a network [l^ . Our 
work reveals that difference of the natural frequency of a 
node with the rest decides its impact on the cluster syn¬ 
chronization and dynamical evolution of all other nodes. 
We present the results for the coupled Rossler oscillators 
on various possible networks such as 1-d lattice, scale-free 
and random networks. Earlier works have shown that the 
network properties, such as degree and betweenness cen¬ 
trality, play an important role in the synchronizability 
of coupled oscillators [l^, and based on the analysis of 
small-world networks it has been shown that the syn¬ 
chronizability can be enhanced by raising the maximum 
degree of the network as well as by reducing the maxi¬ 
mum betweenness [i3- A recent work on power-grid also 
emphasizes on the structural importance of perturbed 
nodes for stability of synchronized state [T^ . 

This Rapid Communication reports that a combina¬ 
tion of the degree and natural frequency mismatch of 
a node with other nodes decides the role of a leader in 
the network, particularly, its impact on the phenomena 
of cluster formation. We demonstrate that, an enhance¬ 
ment in the betweenness centrality does not enhance the 
impact of a leader if the degree is maintained. Apart from 
the impact of a leader on the cluster synchronization, we 
report various different mechanisms of cluster formation 
in the presence of a leader. The earlier works on the 
coupled maps have identified two different phenomena 
for the cluster synchronization, namely, the driven (D) 
and the self-organized (SO) [l^. The SO synchronization 
refers to the state when clusters are formed due to the 
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intra-cluster couplings and D synchronization refers to 
the state when clusters are formed due to the inter-cluster 
couplings. Furthermore, for coupled chaotic oscillators 
having randomly distributed frequencies, it has been re¬ 
ported that with an increase in the coupling strength, 
the nodes with smaller frequency mismatch synchronize 
and form a synchronized cluster while we find that 
if the natural frequency of the nodes are distributed in a 
narrow band, at lower couplings the synchronization be¬ 
tween pair of nodes which are not directly connected is 
preferred. Using a simple star network model, we demon¬ 
strate that the D mechanism is preferred over SO due to 
a common coupling environment D synchronized nodes 
may have. Further, presence of a leader causes an en¬ 
hancement in the synchronization between the directly 
connected nodes thus leading to a transition from the D 
to the mixed or dominant SO clusters state. We study 
phase synchronization of the coupled oscillators, where 
the dynamics of the Uth oscillator can be written as: 


N 


x, = f{x,,i0,) + -J2A^jh{x,,x,y, i = i,...,N (1) 


i=l 


where, Xi S i?"* is m dimensional state vector of the i-th 
oscillator and / : i?™ —>■ i?™ provides the dynamics of 
an isolated oscillator and h is the coupling function. A is 
the adjacency matrix of the network defined as: Aij = 1, 
if oscillators i and j interact, otherwise Aij = 0. The 
degree of a node is defined as ki = P®-" 

rameter e defines the strength of overall coupling among 
the oscillators. 

Coupled Rossler oscillators on networks can be written 
as a set of the following ordinary differential equations: 

Xi = -oJiUi - Zi 

£ ^ 

yi = ujiXi + oj/i -I- — ^ - Vi) 

Zi = b + Zi{xi - c) (2) 

Here uJi is the natural frequency of the oscillator 
which we consider randomly distributed in the interval 
1 < a; < 1.05 0. We take a node acting as a leader 
when its natural frequency is much greater than the rest 
of the nodes in the network (wl :§> 1.05). Later on we 
will explain that the strength of this natural frequency 
mismatch of a node together with the degree of the node 
decide the impact of a leader in a network. In Eq. [21 
a = 0.15, b = 0.4 and c = 8.5 for which the uncoupled 
dynamics is chaotic Q. Further the phase 9 and the 
averaged partial frequency of the z-th oscillator can be 
defined as 6i = arctan^ and fli =< 9i(t) > respectively. 
Here, we note that the frequency Qi is in general different 
from the intrinsic frequency ujj. The uncoupled oscilla¬ 
tors (i.e. e = 0 in Eq. (|T|)) evolve independently. With 
an increase in e, the formation of synchronized clusters 
is observed as e exceeds a critical value £cs- Erequency 
of oscillators may be synchronized forming clusters, i.e. 




FIG. 1: (Color online) (a) and (b) present variation of finter 
(closed circle), fintra (open circle), (fraction of nodes in 
the largest cluster) (open square) and F^ius (fraction of nodes 
forming clusters) (open triangle) with increase in e and ujl 
for the random (a) and scale-free networks (b), (d), (e), (f) of 
N = 50, (fc) = 2. (c) presents the variation of the largest Lya¬ 
punov exponent (A) with change in the coupling strength for 
the random (solid line) and scale-free (dashed line) networks. 
NL, L and PL represent no-leader hub being leader and pe¬ 
ripheral node being leader case respectively. The figures are 
obtained by taking average over 20 random initial conditions. 


U) = U*; j = 1,..., iV;, where is the synchronization 
frequency of cluster I and Ni is the number of oscillators 
in the Tth cluster and I = 1,... ,C; C is the maximum 
number of clusters. 

Depending on the connections between the nodes, rep¬ 
resented by the adjacency matrix, and the synchro¬ 
nized clusters, three phenomena of cluster formation 
have been identified ; driven(D), self-organized(SO), and 
mixed [l^. The quantities fmtra = Nintra/Nc and 
finter = Ninter/Nc, Stand as measures for SO and D 
clusters respectively, where Nintra and Ninter are the 
numbers of intra- and inter- cluster couplings, respec¬ 
tively. In Ninter, Couplings between two isolated nodes 
are not included. Nc is the total number of connections in 
the network. The state which correspond to fintra — 0 
and finter > 0 is defined as the ideal D clusters state; 
fintra > 0 and finter = 0 Correspond to the ideal SO 
state; f^ntra 7^ 0 and finter ^ fintra Correspond to the 
dominant D and finter 7^ 0 and fmtra fmter cor¬ 
respond to the dominant SO clusters state. We take 
fintra — fintra\ < 0.2, to define mixed clusters state 
20|. Eor the higher values, dominant D and dominant 
SO region will shrink and mixed region will grow, while 
for lower values reverse will happen. Eurther, we de¬ 
fine a cluster pattern as a particular state which contains 
information of all the pairs of synchronized nodes dis¬ 
tributed in various clusters in the network. A change in 
the pattern refers to the case when the nodes in the differ¬ 
ent clusters get rearranged [2l[. Eurthermore, we define 
cluster synchronizability in terms of number of the nodes 
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FIG. 2: (Color online) Snap shots showing the change in the 
cluster pattern for random and scale-free networks of = 50, 
(k) — 2 at e = 8. Figs (a), (d) are plotted for random 
network and Figs (c), (d) are plotted for scale-free network. 
The open(red) dots show that the corresponding nodes are 
synchronized and the closed(black) dots show that the corre¬ 
sponding nodes are connected. The snap shot is plotted after 
renumbering the nodes so that the nodes forming a cluster 
come nearby. The star represents the position of the leader in 
the cluster after rearranging the nodes. NL and L represent 
no-leader and leader case respectively. 


participating in a cluster. Based on this, we say cluster 
synchronizability enhances if the number of nodes partic¬ 
ipating in the clusters increases. Further, there might be 
a situation when all the nodes participate in cluster, for 
that cluster synchronizability may be enhanced if the size 
of a cluster increases or the total number of the clusters 
reduces. 

Starting with a set of random initial conditions, we 
evolve the coupled dynamics (Eql2]) on different net¬ 
works, namely, 1-d lattice, scale-free, and random net¬ 
works. After an initial transient we study the cluster 
synchronization. We consider the evolution of coupled 
oscillators without any leader, followed by the investiga¬ 
tion of synchronized clusters in the presence of a leader. 
In the following we discuss the cluster synchronization 
for all the networks. 

The heterogeneity in degree of the scale-free networks 
provides several options for choosing a leader in the 
network yielding very different structural properties to 
the leader. For example, a hub may be assigned as 
a leader making it the highest degree node and conse¬ 
quently best connected with rest of the nodes in the net¬ 
work, or a periphery node may be assigned as a leader 
which makes the leader worse connected with the rest of 
the nodes. As depicted in Fig. [TKa), without a leader the 
coupled dynamics ([2]) leads to the dominant D clusters at 
all the couplings except at very high values where mixed 
clusters exist. For sparse networks considered here, we 
find that while the network exhibits a good cluster syn¬ 
chronization, accompanied with many small clusters, the 
maximum number of nodes in a cluster does not exceed 





FIG. 3: (Color online) Projection of the phase portraits of 
coupled Rossler oscillators for the no-leader and leader case 
for a scale-free network of A = 50, (k) = 2 at e = 2.4. 

20% of the network size. The D clusters correspond 
to the chaotic evolution as reflected by the largest Lya¬ 
punov exponent (Fig. HKc)). What follows that a small 
mismatch in the natural frequencies of the directly con¬ 
nected nodes does not allow them to synchronize with 
each other, whereas the synchronization between a pair 
of nodes which are connected through other nodes gets 
facilitated through a common coupling environment. For 
example, the frequency of the non-identical peripheral 
nodes in a star network synchronize with each other, 
while leaving hub out of the clusters as the dynamics 
of the peripheral nodes from Eq|T]can be written as: 

A, =/(X„cc,)+e(A^-A,); i = 2,...,N 

The hub provides the common coupling to the peripheral 
nodes and thus drive them to form a D synchronized 
cluster. It is however interesting to observe the similar 
behavior for other sparse networks, which consist of many 
star like structures instead of having the ideal situation. 

Upon making a node as a leader by enhancing its natu¬ 
ral frequency higher than that of the other nodes, we find 
that the coupled dynamics leads to a transition from the 
dominant D to the mixed clusters state (Eig. [TJf)). The 
natural frequency of the leader, which leads to this tran¬ 
sition depends on the degree of the node. A hub node be¬ 
ing the leader generates the transition at relatively lesser 
frequency as compared to that required for a peripheral 
node being the leader. For example, for a hub being 
a leader, the mixed cluster state is observed for wl > 2 
(Fig.ITJd)), while for the peripheral node being the leader 
the mixed clusters are observed for wl > 3 (Fig. [TKe)). 

For a hub being the leader, at weak couplings {e < 2.2) 
the clusters remain to be governed by the D mechanism 
(Fig. [TJc)) as observed for the without leader. With an 
increase in the coupling strength, for e: > 3, there is an 
enhancement in the SO synchronization. We emphasize 
that as the D mechanism is still playing a role in the 
cluster formation, with the inclusion of the SO mecha¬ 
nism, the final cluster state becomes of the mixed type. 
Moreover, number of nodes in the largest cluster increases 
(Fig. HKf)). Additionally, in the same coupling regime 
there is a change in the dynamical evolution. The dy¬ 
namics in this regime becomes periodic (Fig. [3Kb) and 
(c)) against the chaotic evolution observed for the no¬ 
leader case (Fig. |3Ka)). 

Another impact of the inclusion of a leader in the net- 
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FIG. 4: (Color online) finter, fintra and fN^ as a function 
of UJL and e. (a), (b) correspond to the random network and 
(c), (d) correspond to the 1-d lattice of (fc) = 4 and N = 50. 
All graphs are plotted for average over 20 realizations of the 
initial conditions. 


work is that it may lead to a completely different clus¬ 
ters pattern (Fig. [2jb) and (d)) than observed for the 
no-leader case (Fig. [JKa) and (c)). 

Inclusion of a peripheral leader at small couplings leave 
the dynamical evolution unchanged. Whereas at strong 
couplings the frequency mismatch of the nodes connected 
with the leader enhances the SO synchronization similar 
as observed for the hub being the leader. This enhance¬ 
ment in the SO synchronization can be explained using 
the revelation that the parameter mismatch between two 
nodes leads to a more stable synchronization [2^ , where 
node having large natural frequency dominates the evolu¬ 
tion with other nodes which are directly connected with 
it leading to the synchronization. 

We remark that the impact of a leader, whose degree 
lie in between the highest and the lowest degree, lies in 
between these two. For example, presence of a leader in 
the ER networks the coupling range 0 < e < 4.0 
keeps the driven phenomena behind the synchronization 
intact, however there is an increase in the SO phenom¬ 
ena. With a further increase in the coupling strength, for 
e > 4.0, the SO synchronization enhances further as in¬ 
dicated by the enhancement in the value of fintra leading 
to the mixed clusters (Fig. SKb)). The enhancement in 
the SO synchronization is associated with the enhance¬ 
ment in the fraction of nodes in the largest cluster as 
depicted in the Fig. HKb). The snap-shots in Fig. [2] de¬ 
pict that inclusion of the leader decreases the number 
of clusters while keep almost all the nodes forming the 
clusters, thereby enhancing the cluster synchronization 
as defined in the model section. 

The 1-d lattice provides an ideal example of a homo¬ 
geneous network. The introduction of a leader in this 
arrangement leads to a similar behavior as exhibited by 
the random and scale-free network depending upon the 



FIG. 5: (Golor online) Snap shots depicting changes in the 
cluster pattern for scale-free networks with N = 50, (k) — 2 
at £ = 3.2. The natural frequency of the nodes are randomly 
distributed in the interval 1 < ca < 1.09 and ojl = 10. 


degree of the lattice which also becomes the degree of 
the leader (Fig. IHd)), which on one hand rules out any 
impact of structural position or preference of the leader 
and on other hand demonstrates importance of the de¬ 
gree of the leader in a network. What follows that the 
impact of a leader on other nodes is high if the leader has 
a large degree, in the absence of which the leader should 
increase the coupling strength in order to bring upon the 
same impact. 

Further, in order to demonstrate the robustness of the 
phenomenon that the introduction of a leader exhibits a 
change in the mechanism of cluster formation, we con¬ 
sider Kuramoto model which is one of the most cele¬ 
brated mathematical model for coupled oscillators [2^. 

N 

6^{t) = ^ X] - ^i)- ( 3 ) 

^ i=l 

We find that introduction of a leader leads to a transition 
from the dominant D to the dominant SO clusters state. 
Fig- El demonstrates formation of the mixed or ideal D 
clusters without a leader and SO cluster in the presence 
of a leader. 

To conclude, we investigate the cluster synchronization 
and phenomena behind the cluster formation for the dif¬ 
fusively coupled Rossler oscillators and find that a leader 
with its natural frequency much higher than that of other 
nodes, has a signihcant impact on the cluster synchro¬ 
nization. The cluster synchronizability of the network is 
enhanced either through an enhancement in the number 
of the nodes participating in the cluster formation or due 
to a merging of several smaller clusters in to larger clus¬ 
ters or due to the formation of larger clusters consisting 
of a completely new set of nodes. Further investigations 
reveal that the introduction of a leader may lead to a 
transition from the D to SO mechanism of the cluster 
formation. Thus, in the presence of a leader, synchroniza¬ 
tion between the directly connected nodes is enhanced. 
The presence of a leader may also lead to a transition 
from the chaotic to the periodic dynamical evolution and 
hence a leader may be introduced for chaos control iby 
tuning frequency of a single node. Interestingly, a leader 
has the maximum impact on the cluster synchronization 
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if it is the highest degree node in the network rather 
than being the node which has the highest betweenness 
centrality. If a leader has a lower degree, its natural fre¬ 
quency should be relatively higher in order to achieve 
the enhancement in the cluster synchronization. For ho¬ 
mogeneous networks, where all the nodes have the same 
degree, coupling should be high in order to have a tran¬ 
sition from the dominant D to the mixed cluster state. 
Furthermore, the presence of a leader not only changes 
the phenomena behind the cluster formation but may 
also completely change the cluster pattern. 

Leaders naturally arise in real-world networks such 
as in social networks [l^, neural networks [l^, protein 
translation regulatory networks [1^. In social network, 
a leader may possess one of the characteristics of power. 


experience, fame or wealth, while in biological networks, 
such as in neural and the protein translation regulatory 
networks, a leader may be one which initiate certain pro¬ 
cesses [ij, [15 . Our work may be extended to capture 
particular properties of a leader for understanding the 
origin of synchronized clusters in these systems [ig. For 
example, a leader may have different coupling strength, 
such as in the brain network, where the synapses become 
weak with age Big. 
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